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We determine the improvement factor c$w in one-loop lattice perturbation theory for the plaquette 
and Symanzik improved gauge actions. The fermionic action is 0(a) clover improved with one- 
time stout smearing, csw is derived from the one-loop correction to the quark-quark-gluon vertex 
in the off-shell regime. We give a first numerical value for the one-loop contribution to the non 
gauge-invariant improvement coefficient c^gi f° r the quark field using the plaquette action. A 
discussion of mean field improvement is included. 
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1. Introduction 

Current simulations with 2 + 1 flavors require highly improved gauge and quark actions. 
Renormalization group improved gauge actions are order £?(a 2 ) improved and should be preferred 
to the @{a) improved plaquette gauge action. In accordance with our numerical simulations we 
take the Symanzik improved gauge action [p], §1 

s Syman Z ik = ^ [ C() i Re Tr ( 1 - U plaquette ) + d £ rm ± ReTr (1 - U rectangle )] (1.1) 

with ^ 

d = - — , c = l-8ci. 

As the fermionic action, we use the clover improved action as proposed by Sheikholeslami and 
Wohlert [Q] which means that one has to add the so-called clover term to the standard Wilson 
fermion action 

Sf' ver = Sf lson - c sw £ £ ig r -% CT MV F M v(«) Wn , (1-2) 

n H,V 

where F^ v (n) is the field strength in clover form and o^y = ij2{y^y v — y v %i)- An additional im- 
provement can be achieved with ultraviolet filtering or smearing the gauge links Uu in the fermionic 
Wilson action Sp Son : it reduces the chiral symmetry breaking of Wilson quarks among light fla- 
vors. There have been proposed several smearing techniques (for a detailed discussion see [Q]). We 
use the stout smearing of Morningstar and Peardon It is given by a sequence of transformations 

U^U^ -^u£ ) ---^U { ; ) =%, (1.3) 

with 

U { ; l+l) (x)=e i Q i ^ u -^U^(x). 

The function <2u (U, (»^ v ) depends on the staples of the gauge link under consideration and on the 
stout parameters (0^ v which determine the strength of smearing. We chose an isotropic parameter 
ft)^v = CO and one step smearing which is recommended by various investigations. 



It is of importance to determine the improvement factor c$w appearing in ( |1.2| ) as precisely 
as possible. Non-perturbative determinations are to be preferred but for the combination described 
above there are no results obtained so far. In perturbation theory c$w has the form 

c sw = l+g 2 c ( s ^ + 0(g 4 ). (1.4) 

There have been published results for for plaquette action with twisted antiperiodic boundary 
conditions |Q] and Schrodinger functional method [j^]. For some popular improved gauge actions 
Aoki and Kuramashi [||] calculated the one-loop correction using conventional perturbation theory. 
All results are obtained for unsmeared gauge links in the on-shell regime. 

In this paper we calculate for Symanzik improved gauge action with stout smearing in 
conventional perturbation theory. We do the calculation off-shell. This enables us to determine 
the one-loop contribution to the non gauge-invariant improvement coefficient c^gi for the quark 
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fields y as proposed in Using BRST symmetry arguments the authors proposed the off-shell 
improvement for the quark fields 1/4 to be 

Y* = (l+ac D + aigc NGI 4)Y, (1-5) 

where the coefficient cngi does not contribute on-shell. Its perturbative expansion is known to 
be Q] 

CNGI=g 2 C { K GI + 0(g 4 ). (1.6) 

In order to determine cjj!^ either a two-loop calculation of the quark propagator or a one-loop 
calculation of the quark-quark-gluon vertex is required.. The improvement coefficient cd has been 
calculated to one-loop order in [[R]] . 

2. Improvement procedure 

In the approach of conventional perturbation theory we use the quark-quark-gluon vertex 
A^(pi,p2) as discussed in [Q] already. Looking at the ff(a) expansion of tree-level i{pi,pz) as 



derived from action (1.2) 

^(P2,Pi) = -igj^-gjarl(pi +P2)ii-c S wigjara^ a {P2-P\)a + ^{a 2 ), (2.1) 



one can see by inserting (1.4) that a one-loop calculation for An(p2,pi) provides necessary condi- 



tions to determine c^. We omit in all three-point functions the common overall color factor t ac . 
In ( |Q| ) p\ (P2) are the incoming (outgoing) momenta. The off-shell improvement condition states 
that the non-amputated improved three-point function G* yil {p2,p\) has to be free of 0(a) terms in 
one-loop. The unimproved and improved three-point functions are defined by 

Gn(P2,Pi) = S(pz)Ay(p2,pi)S{pi)D Vfl (q), (2.2) 
G* >fl (j>2,pi) = S*(p 2 )A*,v(P2,pi)S*(pi)D V n(q), (2.3) 

with q = P2 — Pi- D v ^(q) is the full gluon propagator which is G (a) -improved already. A^(p2,pi) 
and A*n(p2,pi) are the unimproved and improved amputated three-point functions. The corre- 
sponding quark propagators are given by 

S- l (p) = i^ P {p) + ^wip) = i^ P {p)U-\ai]^^^\ , (2.4) 
S- l (p) = ijl. p {p). (2.5) 



In terms of the improved quark fields (1.5) G^(p2,pi) can be related to its improved version 



G^p2,Pi) = G^(p2,Pi)-aigc NGI ^[((4^- l + ]^- 1 4) Ap)] . (2.6) 



In deriving (2.6) we have assumed (A) = 0, & denotes the Fourier transform. Taking into account 



l.q ) we insert in our one-loop calculation the corresponding tree-level expressions 

aigc NGI &[({4l?- l +y- x 4)A») tree ] =*ig 3 4 G i(yv^r + ^r7 v ) D% e (q), (2.7) 

\ 1 Pi 1 F2 / 
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or its amputated version 



aigc NG] & [({tfp-' + I?- 1 4) A M >££] = -ag 3 c% (/ 2 ft + Jn A 



(2.8) 



If we amputate (2J2) and use ( |2.4| ), (2.6) and ( |2.8| ) we get the off-shell improvement condition 



A / i(p2,/>i) = A* )M (p 2 ,Pi) + ag 



3 Ji) 



c 



^7m + 7m^ 



2"-^£ p ( P2 

which should hold to order '■ 



by determining c^ 7 and c^J, correctly. 



(2.9) 



3. Calculation 



q=P 2 -P 1 






(d) 





Figure 1: One-loop diagrams contributing to the amputated quark-quark-gluon vertex 

The diagrams contributing to the amputated one-loop three-point function are shown in Fig. []]. 
The calculation is performed combining symbolic and numerical methods. For the symbolic com- 
putation we use a Mathematica package that we developed for one-loop calculations in lattice 
perturbation theory (for a more detailed description see [JTTJ]). It is based on the infinite volume 



algorithm of Kawai et al. ||10p. The analytic treatment has several advantages: one can extract the 
infrared singularities exactly and the results are given as functions of lattice integrals which can 
be determined with high precision. The disadvantage consists in very large expressions especially 
for the problem under consideration. In the analytic method the divergencies are isolated by dif- 
ferentiation with respect to external momenta. As can be seen in Fig. [j] diagrams (b) and (c) have 
two gluon propagators. So no parametrization can be chosen with only internal momentum flowing 
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through the gluon lines. Therefore at least one gluon propagator has to be differentiated. Look- 



ing at the full analytic form of the gluon propagator for improved gauge actions Q13| ] one easily 



sees that huge analytic expression would arise. As discussed in [ |13| ] one can split the full gluon 
propagator D n ^ mved (k) 

D taproved (jt) =D plaquette (jt)+AD ^ (jk) (3J) 

The diagrams with v quette (k) only contain the logarithmic parts and are treated with the analytic 
Mathematica package. The diagrams with at least one AD jlv (k) are infrared finite and can be de- 
termined safely with pure numeric methods. We have written a C program with a Gauss-Legendre 



integration algorithm in four dimensions (for a description of the method see [ ]1 1 , |12[]). We choose 
a sequence of small external momenta {p\,pi) and perform an extrapolation to vanishing momenta 
in order to extract the corresponding values. Additionally, we have written an independent FOR- 
TRAN code which computes the one-loop contributions for each diagram including the infrared 
logarithms. Results for both methods agree within accuracy. 

The Feynman rules for non-smeared Symanzik gauge action have been summarized in [§]. For 
the stout smeared gauge links in the clover action the rules are given for the forward case by [Q]. The 
corresponding Feynman rules needed for the quark-quark-gluon vertex are much more complicated 



and have been derived by the authors. They are too long as to be given in this proceedings H15Q. 

The calculation has been done in Feynman gauge with Wilson parameter r = 1. All the one- 
loop coefficients are calculated at csw = 1 because g 3 csw = g 3 + 0(g 5 ). 

4. Results 

The anticipated general structure for the amputated three-point function in one-loop is 

A M (p2,Pi) = tf[ S {P2,Pi)+Aiatig 3 7n 

+Biat ^g 3 (jf 2 7n + 7n P'x) + Qat y g 3 <y^ a q a (4. 1) 

A^f {p2,pi) is the universal part of the three-point function independent of the chosen gauge action 
computed in the MS-scheme 

a , x a 



A M (j)2,pi) = -igYn-g^l{phfi+P2,n) -c sw ig^ a q a 

+ig 3 Fi,n(pi,P2,q)+a g 3 F 2 ^ (pi,p 2 ,q). (4.2) 



F\,ii{p\,pi,q) and F2u(pi,P2,q) are complicated functions involving poly logarithms and loga- 
rithms. They will be given in [O]. The quantitites Ai at , B\ at and Ci at are obtained as 



A, at = C F (0.03783- 0.93653 (0 + 3.42833 co 2 + 0.01266 log(aji)) 

+N C (-0.02200 + 0.01266 log(a/x)) , 
Bi a t = C F (0.03804 - 1. 03749 co + 3.43791 co 2 +0.02533 log(aju)) 

+N C (-0.02432 + 0.01925 co + 0.01266 log(aju)) , (4.3) 
Ciat = C F (0.11618 + 0.82813(0 -2.45508 co 2 ) 

+N C (0.01215 + 0.01109 w -0.30228 co 2 ) , 
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with Cp = (iV~ — 1 ) / (2N C ) for SU (N c ). As shown in ( 2.S ) we need the self energy parts L p (p) and 



L\v (p) as defined in ( |2.4j ) to solve the off-shell improvement condition 



16k 2 

*w(p) = 1 " fg|r [2 log(ap) 2 +L 2 ] . (4.4) 



i » 



It turns out that the self energy parts Ei and £2 contribute only to c^ GI . For the Symanzik gauge 
action we will present them in [|i5|]. For the plaquette action we get 

L plaq = 8 20627 - 196.44600 co +739.68364 co 2 , 

l2! aq = 7.35794-208.58321w+711.56526w 2 . (4.5) 



We use ( |4.2| ) and d4.4[ ) to construct the left hand side of ( |2.9[ ) whereas ( |4.4| ) with ( |4.5[ ) are inserted 
into the right hand side. In order to fulfill we get the following improvement coefficients for 
the plaquette action 

c {\,piaq) = Nc (0.00143 - 0.01 166 <o) , (4.6) 
c {^piaq) = Cp (0.16764+ 1.07915 w-3.686680) 2 ) 

+N C (0.01502 + 0.00962 co - 0.28479 co 2 ) . (4.7) 

For the Symanzik improved gauge action we find the improvement coefficient 

= c F (0.11618 + 0.82813W -2.45508m 2 ) 

+N C (0.01215 + 0.01 109 co- 0.30228 co 2 ) . (4.8) 

5. Mean field improvement 

It is known that lattice artefacts make the perturbative expansion worse. One possible im- 
provement procedure is to replace the naive coupling constant g by its mean field improved value 
gMF = g/uf) where is the average plaquette value for the corresponding gauge field action. By 
scaling all gauge links in the clover field strength F^ v (n) in ( |1.2[ ) by 1/uq one obtains the mean 
field improved csw as 

c jw="o c w- (5.1) 
The perturbative expansion of uq is known to be 

«o = l-%^ H , (5-2) 

where k u for popular gauge actions are given in [O]. Therefore, the perturbative expression for the 
mean field improved csw is given by 

csw = c^u^ = ^ (l+dfF^-^^+^Cglo,)) =<%Z* + 0<&f)- (5-3) 
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For the future simulations of the QCDSF collaboration we have the following numbers for the 
Symanzik action and 2+1 flavors 

C F =4/3, N c = 3, i4 = 0.6065, gjf F = 1.71335, k u = 0.732524 n 2 . 
This gives the one-loop expression for csw parameter as 

c sw = l+g 2 (0.19136 + 1.13745w-4.18029co 2 ) + ^(g 4 ), (5.4) 
c sw P = 4 ( l + Smf (0. 19136 + 1. 13745 co - 4.18029 co 2 )-g 2 MF 0.183 13) 

= 1.47557 + 2.83568 w- 10.42148 co 2 (5.5) 



For no stout-smearing (co = 0) the result ( p.4| ) has to be compared with the number given in [|8J] : 
C SW^ = 0.19624449(1). The minor difference to our value c™ = 0.19136 can possibly be related 
to an inaccuracy in our numerical integrations. In the simluation the stout parameter co is chosen 
to be co = 0.1 leading to a mean field improved value c^ p = 1.65492. 



This investigation has been supported by the DFG under contract FOR 465 (Forschergruppe 
Gitter-Hadronen-Phanomenologie). 
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